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ABSTRACT 

We observe a solar surge in NOAA AR11271 using SDO AIA 304 A image 
data on 25 August, 2011. The surge rises vertically from its origin upto a height 
of :^65 Mm with a terminal velocity of ^100 km s~^, and thereafter falls and 
fades gradually. The total life time of the surge was ^20 min. We also mea- 
sure the temperature and density distribution of the observed surge during its 
maximum rise, and found an average temperature and density of 2.0 MK and 
4.1x10^ cm~^, respectively. The temperature map shows the expansion and 
mixing of cool plasma lagging behind the hot coronal plasma along the surge. 
As SDO/HMI temporal image data does not show any detectable evidence of 
the significant photospheric magnetic field cancellation for the formation of the 
observed surge, we infer that it is probably driven by magnetic reconnection 
generated thermal energy in the lower chromosphere. The radiance (thus mass 
density) oscillations near the base of the surge are also evident, which may be 
the most likely signature of its formation by a reconnection generated pulse. In 
support of the present observational base-line of the triggering of the surge due 
to chromospheric heating, we devise a numerical model with conceivable imple- 
mentation of VAL-C atmosphere and a thermal pulse as an initial trigger. We 
find that the pulse steepens into a slow shock at higher altitudes that triggers 
plasma perturbations exhibiting the observed features of the surge, e.g., terminal 
velocity, height, width, life-time, and heated fine structures near its base. 

Subject headings: magnetohydrodynamics (MHD) : sun — corona: sun — chromosphere 
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Introduction 

Various types of plasma ejections are significant in t he solar atmosphere to transpo rt 



mass and enerev at diverse soatio-temooral scales fe.e.. iBohlin et al. 


1975; 


Innes et al. 


1997; 


Shibata et al. 


2007 




Katsukawa et al. 


2007; 


Culhane et al. 


2007 


FilioDOv et al. 




2009; 


Srivastava & Murawski 


2011; 


Morton et al. 


2012. and references cited there). The 



magnetically structured lower solar atmosphere is dynamically filled with diflFerent types 
of jet-like phen omena at short spatial scales, e.R., spicules, mottles, fibrils, anemone jets, 



tornadoes (e.g.. 



De Pontieu et 



Wedemeyer-Bohm et al.l 



al. 



2004 



Hansteen et al. 



200^; 



De Pontieu et al. 



20074 



2011 



2012al . and references cited there). The large-scale plasma jets 



(polar jets, s urges, sprays etc) are a. 



corona (e.g., 



Srivastava fc Murawski 



so important to channel 



2011 



Uddin et al. 



2012 



mass and energy into upper 
, and references cited there). 



The magnetic reconnection and magnetohydrodynamic (MHD) wave activity are found 



Yokovama fc Shibata 



to be the basic mechanisnis for driving such plasma eject a in the so. 



1995 



19961: 



Innes et al...l997.: .De Pontieu et al 



ar atmosphere (e.g.. 



20041: 



Hansteen et al. 



2iM 



Shibata et al.l 120071: iKatsukawa et al.l 120071: 



Cirtain et al.l l2007l: iDe Pontieu et al. 



Kamioetal 



20101 : 



2007alJblJcl: iPariat et al.l l2010l: iFihppov et al.l l2009l: iMuraws 



Murawski et al 



2011 



Morton et al. 



d fc Zaqarashvilil l2010l : 



20121, and references cited there). 



Among various types of solar jets, the solar surges are cool jets, typically fo rmed by a 



plasm a that is usually visible in H^^ and other chromospheric and coronal lines ( ISterling 



20001 ) . The solar surges are most likely triggered 



by magnetic recon nection and they 



20071 ). Solar surges are 



exhibit diflFerent episodes of heating and cooling ( iBrooks et al. 
mostly associated with fiaring regions and places of transient and dynamical activities in 
the solar atmosphere. The solar surges and their evolution at various temperatures have 



been observed extensiv ely in association with magnetic field and fiaring activities of t 



solar atmosphere (e.g.. 



Schmieder et al. 



1994; 



Chae et al. 



im 



ne 



Yoshimura et al. 



200^ 
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Liu fc Kurokawal l2004l : Ijiang et all 120071 : lUddin et all 120121 . and references cited there) 



Inspite of th e direct magnetic reconnection and photospheric 



cancellation ( IGaizauskas 



1996 



Sterling 



20001 : 



Liu et al. 



20051 : 



magnetic fie 



Uddin et al. 



surg:es can also be t riggered by impulsive generation of pressure pulse (IShibata et al 



d eme rgence and 



2012h. the solar 



Sterhng et al. 



19931 ). and may also be associated with the explosive events ( IMadjarska et al. 



1982 



20091 ). However, understanding of their exact drivers requires further investigations, and 



more than one mechanism may work for the surge evolution depending upon the local 
plasma and magnetic field conditions. Therefore, study of the driving mechanisms of any 
particular surge is always the forefront of solar research that may also be useful probing of 
the localized conditions and magnetic activities in the solar atmosphere. 

In the present paper, we describe empirical results derived from the observations of a 
surge that originated from the western boundary of an active region NOAA AR11271 on 
25 August, 2011. We study the kinematics and plasma properties of this surge, as well as 
its most likely trigger mechanism. We select a distinct surge that was not related with the 
recurrent surge activities in AR11271, and apply a numerical model that reproduces its 
observed properties. Our model suggests that the trigger of the observed surge may be due 
to the reconnection generated thermal pulse in the lower solar atmosphere. In Sect. 2 we 
describe the observations of the surge, its kinematics, and derived plasma properties. We 
also investigate the magnetic fields at the base of the surge using Helioseismic Magnetic 
Imager (HMI) data and discuss its consequences. We report on the numerical model of the 
observed surge and present numerical results in Sect. 3 and conclude with the discussions. 
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2. STUDY OF THE KINEMATICS, PLASMA, AND MAGNETIC 
PROPERTIES OF A SOLAR SURGE TRIGGERED FROM AR11271 

For this study we have used the imaging and magnetic field o bservations fr om 



2OI2I) 



instruments onboard SDO. AIA is a multi-channel full-disk imager ( iLemen et al. 
with resolution elements of 0.6" and a cadence of 12 s. For this study we have used data 
of 171 A (coronal) and 304 A (upper chromospheric/TR) filters. HMI is designed to 
study the photospheric magnetic field and takes full-disk images at a cadence of 3.75 s to 
pr ovide Doppler inte nsity, magnetic field, and full vector magnetic field at a cadence of 45 



s (Schouetal 



20121). SDO/AIA and SDO/HMI data have been obtained from Solar-Soft 
cut-out service established at Lokheed Martin Space Research Laboratory, USA Q, which 
are corrected for fiat-field and spikes. For further cleaning and calibration of this data, we 
also run aia_prep subroutine of SSWIDL. 

We use an image sequence of a solar surge in AR11271 on 25 August 2011 during 
09:11-09:31 UT as observed in the 304 A and 193 A filters, which represent 10^ K and 
10^ K plasma, respectively. The surge originated near the western boundary of AR11271 
as shown in Fig. 1. These images are overlaid by SDO/HMI magnetic field contours at a 
level of ±800 G : yellow (blue) stays for positive (negative) polarity.lt is clearly evident 
that the surge occurred near a positive polarity region (yellow) in the north-west side of 
a big negative polarity spot (shown by blue circle). The surge plasma mostly emits in 
the 304 A filter sensitive to the plasma temperature around ^10^ K , but does not show 
much emission in the coronal filter. However, some emission in the hotter coronal channel 
is visible up to the height of 20 Mm from the foot point of the surge (cf., right panel of 
Fig. 1). Therefore, the surge consists of multi-temperature plasma. 



ht t p : / / www . Imsal . com / get _aia_dat a / 
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2.1. Kinematics of the Surge 

Fig. 2 represents the spatio-temporal evolution of the surge material using the sub-field 
data. It is clear that the surge is being activated at 09:11 UT near the western boundary 
of NOAA AR11271, and faded around 09:21 UT. Th e observed sur^e was not associated 



2OI2I ). This image 



with the recurrent spray surge activity in AR11271 ( iKayshap et al. 
sequence shows that the surge reaches upto a projected height of ^65 Mm in ?^660 s. The 
terminal speed of the surge is found to be ?^100 km s~^. This surge is well collimated with 
a maximum width of ^7 Mm and does not show any significant lateral expansion. Its 
down falling time is ^^500 s in which it also faded gradually. Top-panel of Fig. 3 shows 
the distance-time measurement along the observed surge. It is evident that the denser 
core material of the surge first rises up and then decelerates. Some fainter and less denser 
plasma reaches much higher into the solar atmosphere, and falls freely thereafter. The 
distance-time map across the surge at a particular height around ^15 Mm is shown in the 
bottom panel. The radiance periodicity is observed on 3.0 min time-scale, which may be the 
signature of the periodic enhancement of mass density in the surge. This may provide the 
clues about the periodic driver of the surge, which will be discussed in detail in Sec. 3. It is 



also possible that the periodicity of 3.0 min near the sur ge base may be due to the swir 



motion as recent 



De Pontieu et al. 



y obs erved in spectroscopic data (e.g.. 



Wedemeyer-Bohm et al 



2012b 



mg 



20121, and references cited therein). However, such small-scale motion is 
below the detection limit in our data. Moreover, the radiance oscillations occur higher in 
the transition region, and therefore, in the present case it may be the most likely evidence 
of the periodic signature of the arrival of slow shocks in the upper chromosphere/transition 



region. 
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2.2. Estimations of the Average Temperature and Density of the Surge 



In the present section, we carried out emission measure and temperature distribution 
analysis, as well as the tracing and estima tion of density aloi i R the sur^e by extensively 



using the automated method developed by 



Aschwanden et al. 



(1201 ll ). Understanding such 



plasma properties provide more clues about the driver of the surge, while comparing with 
the numerical results. SDO/AIA data provides infor mation of the emission m easure and 



estimation of the average density and temperature (Aschwanden et al, 



201lh . We use 



full-disk SDO/AIA in all EUV channels at the time of the maximum rise of the surge 
around 09:22 UT on 25 August, 2011. We calibrate and clean the data using aia_prep 
subroutine of SSW IDL as well as co-align the AIA images as observed in va rious AIA 



filters by using the co-alignment test as described by 



Aschwanden et al. 



(1201lh . 



Using the automated code developed by 



Aschwanden et al. 



( 120111 ). we obtain emission 



measure and temperature maps for six AIA filter full-disk and co-aligned images (304 A, 
171 A, 193 A, 94 A, 335 A, 211 A) in the temperature between 0.5-9.0 MK. The emission 
measure map in the top-left panel of Fig. 4 shows that the surge has significant emissions 
from its denser core part. The emissions from its leading edge is comparatively low. 
Top-right panel of Fig. 4 shows the distribution of temperature in and around the surge 
location. It is evident that heating is confined near its foot-point. Greenish- yellow regions 
represent the highest temperature in and around the footpoint of the surge. The reddish 
region shows the plasma at the base of the surge with a temperature of 6.2-6.3 MK. Above 
the heated area near the base of the surge, the cool plasma is maintained at sub-coronal 
temperatures as represented by blue-black colors. It is also mixed with the comparatively 
hot plasma. The leading edge of the surge is, however, maintained at comparatively high 
coronal temperatures. The heating near the base of the surge probably drives it, and the 
cool plasma moves behind the hot plasma at its leading edge. It is also mixed with the hot 
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plasma, and constitutes the multi-temperature surge. 



Using the automated procedure (lAschwanden et al. 



20111), we trace the loop segments 



and surge in the co-spatial and co-temporal images of various wavelengths of SDO/AIA 
data. The input free parameters for this automated tracing of the surge is summarized in 
Table 1. The bottom panel of Fig. 4 shows that the surge material is erupting along open 
field lines, indicated by white arrow, from the western part of active region boundary. The 
surge consists of multi-temperature plasma. The contribution of cool surge plasma (T=0.1 
MK) as evident in 304 A channel, is represented by yellow colour. For the estimation of 
the plasma properties (temperature and density) along the surge material and its overlying 
diffuse corona, we trace the larger segment of this coronal structure using another set of free 
parameters (cf.. Table 1) and implement the best forward fit of the possible DEM solutions. 
We choose this set of free parameters on the basis of trial and error, and it gives the best 
DEM forward fit in the case of 171 A fiuxes at each pixel of the considered region around 
the surge and its overlying atmosphere. In Fig. 5, the estimated density and temperature 
are displayed. The chi-square of the estimation in the forward modeling is good-enough 
up to 70 Mm along the surge as well as upto 140 Mm in its overlying diffuse corona. The 
average temperature and density of the surge material are respectively 2.0 MK and 4.1x10^ 
cm~^. We can estimate these physical parameters in all the traced structures in all AIA 
channels. However, we display the best fit results of the estimation of the density and 
temperature along the surge and its overlying corona in 171 A channel as shown in Fig. 5. 



2.3. SDO/HMI Investigations of the Magnetic Field Properties at the Base of 

the Surge 



The surge occurred at the north-western side of the big sunspot with negative polarity. 
It was located at the outer layer of a positive polarity region (cf., snapshot on 09:00 UT 
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in Fig. 6). The surge is shown as intensity contour (blue) of co-ahgned AIA data on HMI 
snapshot of 09:21 UT in Fig. 6. The morphol ogy of the or i gin si te of the observed surge 
is almost similar to that observed recently by 



Uddin et al. 



( 120 12l ) in association with the 



recurrent surges from NOAA 10884 on 25 October 2003. However, we do not observe any 
strong signature of the magnetic field annihilation as well as moving magnetic fe a tures 



(MMFs) in our present observational data as previously observed by .Uddin et al. 



fl2012h . 



Therefore, we infer that the observed surge may not be triggered by strong photospheric 
magnetic field activities. It is evident that tiny negative magnetic polarities have been 
evolved near the surge location, as well as the major positive polarity region has shown some 
areal enhancement towards the western direction. Therefore, such rearrangements of the 
magnetic polarities a round the surge region may pro duce more complexity in the overlying 



magnetic fiux-tubes (IRoy 



1973: 



Canfield et al. 



19961 ). e.g., in low- lying loops and an arch 



crossing near the base of the surge (cf., 09:11 UT AIA 304 A snapshot of Fig. 2; HMI- AIA 
snapshot on 09:20 UT in Fig. 6). Therefore, the chromospheric reconnection between the 
open field lines as well as surrounding closed fields may trigger the surge. The localized 
brightening and heating at the same location is evident where dark arch type fiux-tube and 
low- lying closed loops are crossing near the surge base (cf., 09:11 UT snapshot of Fig. 2, 
and right panel of Fig. 4), which may generate the thermal pulse probably triggering the 
observed plasma perturbations. In the next section, in support of the present observational 
base-line of the triggering of the surge due to chromospheric heating, we devise a numerical 
model with implementation of VAL-C atmosphere and a thermal pulse as an initial trigger, 
which re-produce the observations closely. 
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3. A NUMERICAL MODEL FOR THE PULSE-DRIVEN SOLAR SURGE 

To reproduce the observed surge, our model system assumes a gravitationally-stratified 
solar atmosphere that is described by the ideal two-dimensional (2D) MHD equations: 

^ + V-(^V) = 0, (1) 



Q^ + g{y-V)\^-Vp+^{VxB)xB + gg, (2) 



^ + V-(pV) = (l-7)pV-V, (3) 



_ = Vx(VxB), V.B = 0. (4) 

Here g, V, B, p = ^^T, T, 7 = 5/3, g = (0, -g) with its value ^ = 274 m 8"^ m, feg, 
are respectively the mass density, flow velocity, magnetic fleld, gas pressure, temperature, 
adiabatic index, gravitational acceleration, mean particle mass, and Boltzmann's constant. 
It should be noted that we do not invoke the radiative cooling and thermal conduction 
in our model for simplicity reasons since we simulate the cool surge eject a in the solar 
atmosphere. 

3.1. Equilibrium configuration 

We assume that at its equilibrium the solar atmosphere is still (Vg = 0) with a 
force-free magnetic field, 

(V X Be) X Be = , (5) 



such that it satisfies a current-free condition, V x Bg = 0, and it is specified by the magnetic 
flux function, A, as 

Be = V X (Az) . (6) 

Here the subscript 'e' corresponds to equihbrium quantities. We set an arcade magnetic 
fleld by choosing 

A{x, y) = BqAb cos {x/AB)exp[-{y - 2/r)/AB] (7) 
with Bq being the magnetic field at ^ = ^r, and the magnetic scale-height is 

2L 



(8) 



We set and hold fixed L = 100 Mm. For such choice the magnetic field vectors are weakly 
curved (Fig. [71 bottom panel). 

As a result of Eq. (Ej) the pressure gradient is balanced by the gravity force. 



(9) 



With the ideal gas law and the t/-component of Eq. 1^, we arrive at 



Peiy) =Pq exp 



^ dy 



Qe{y) 



Pejy) 
gHy) 



where 



My) = 



mg 



(10) 



(11) 



is the pressure scale-height, and po denotes the gas pressure at the reference level that we 
choose in the solar corona at = 10 Mm. 



We take an equilibrium temperat ure profile Tp(^) 



from the VAL-C atmospheric model of 



Vernazza et al 



'or th e solar atmosphere derived 



(jl98ll ). Having specified T^^y) (see 



Fig. 7, top panel) with Eq. ( ITOl) we obtain the corresponding gas pressure and mass density 
profiles. 
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The transition region (TR) is located y 2.7 Mm. Above the TR we assume an 
extended solar corona. Below the solar chromosphere, the temperature minimum is located 
dityc^ 0.9 Mm. 

3.2. Pert urbat ions 

We initially perturb the equilibrium impulsively by a Gaussian pulse in a gas pressure 



VIZ., 



p{x,y,t = 0) =Pe{y)+Pe{y)A^exp 



(x - xo)^ {y - yof 



(12) 



Here is the amplitude of the pulse, (xo,?/o) is its initial position and Wy^^ Wy denote its 
widths along the x- and ^/-directions, respectively. We set and hold fixed = 80, Xq = 
Mm, yo = 1.75 Mm, = 2.25 Mm, and Wy = 0.2 Mm. 

3.3. Results of the Numerical Simulations 



Equations ffl -dH) are solved numerically using the FLASH code ( ILee fc Deand 12009 



Fryxell et al. 



2000l ) . This code uses a second-order unsplit Godunov solver with various 



slope limiters and Riemann solvers, as well as adaptive mesh refinement (AMR). We set the 
simulation box of (—12.5, 12.5) x (0, 125) Mm^ along the x- and ^-directions (Fig. El bottom 
panel) and set and hold fixed all plasma quantities at all boundaries of the simulation region 
to their equilibrium values which are given by Eqs. (|6j) and f fTOl) . These fixed boundary 
conditions performed much better than transparent boundaries, leading only to negligibly 
small numerical refiections of wave signals from these boundaries. 

In our modeling, we use an AMR grid with a minimum (maximum) level of refinement 
set to 4 (7) (FigE]). The refinement strategy is based on controlling numerical errors in 
mass density, which results in an excellent resolution of steep spatial profiles and greatly 
reduces numerical diflFusion at these locations. 



Fig. 9 displays the spatial profiles of the plasma temperature (color maps) and velocity 
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(arrows) resulting from the initial pulse of Eq. ( fT2l) . which splits into counter-propagating 
pulse trains. The downward propagating part is reflected from the dense ambient plasma 
of the photosphere. This reflected part lags behind the originally upward propagating 
pulses, which become slow shocks. Because the plasma is initially pushed upward, the 
under-pressure results in the region below the initial pulse. This under-pressure sucks up 
comparatively cold chromospheric plasma, which lags behind the shock front at higher 
temperature. As a result, the pressure gradient force works against gravity and forces the 
chromospheric material to arrive into the solar corona in form of a rarefaction wave. At 
t = 200 s this shock reaches the altitude oi y 65 Mm and the cold plasma edge is located 
Sit y 35 Mm, which is below the slow shock. The reason for the material being lifted up is 
the rarefaction of the plasma behind the shock front, which leads to low pressure there. The 
next snapshot (top middle panel) is drawn for t = 400 s. At this time the chromospheric 
plasma is located below Sit y 55 Mm. At the next moment, t = 600 s (top right panel) the 
chromospheric plasma shows its developed phase, reaching the level of ^ = 65 Mm, which 
approximately matches the observational data of Fig. 2 (top right panel). The cool surge 
slows down while propagating upward. At t = 800 s it arrives at the level of y 70 Mm 
and subsequently the surge already subsides being attracted by gravity (bottom panels). 
At the subsiding stage the heated flne structures develop as a result of interaction between 
the down-falling plasma. During the evolution of the surge, the comparatively hot plasma 
is evident near the base of the simulated surge that qualitatively matches the observations. 
The leading front of the slow shock moving above the cool surge plasma may also leave 
some traces in the overlying diflFused quiet-Sun atmosphere evident in the coronal fllters (cf., 
right panel of Fig.l). However, the less denser and diffused overlying quiet-Sun may not be 
visible as a tracer of a slow shock front. The temporal evolution of the density reveals that 
the surge is denser near its base, while the density decreases with the height (Fig. 10). This 
is consistent with the observational data. 
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4. DISCUSSIONS AND CONCLUSION 

In the present paper, we study the kinematics and plasma properties of the solar surge 
triggered in AR11271 on 25 August, 2011. We also performed a 2 D numerical simulation 
usinR VAL-C mode l of the solar atmosphere as initial condition and the FLASH code 



( IFryxell et al 



20001 ) . which approximately reproduce the dynamics of the observed cool 
surge. The surge rises vertically from its origin site up to a height ?^65 Mm with a terminal 
velocity of ?^100 km s~^, and thereafter falls and fades gradually within its total life time of 
^20 min. We also measure the temperature and density distribution of the observed surge 
during its maximum rise, and its average temperature and density are respectively 2.0 MK 
and 4.1x10^ cm~^. We do not get any strong evidence of the photospheric magnetic field 
activity (i.e., MMFs, fiux-emergence, field annihilation, etc) in support of the formation of 
the observed surge. The surge is most likely found to be driven by magnetic reconnect ion 
generated thermal energy in the chromosphere. 

The numerical simulations exhibit a general scenario of rising and subsequently 
down-falling/fading of the observed surge with a number of observational features such as 
its height and width, ejection, approximate rising time-scale, and average rising velocity 
that are close to the observational data. The heated fine structures at the base of the 
surge are also matching well both in the observations and theory. The temperature 
distribution (Fig. 4) shows that the significant amount of cool plasma expands along the 
surge, while its leading edge is maintained at coronal temperatures. This scenario matches 
qualitatively with our simulation results when the comparatively cooler core plasma of 
surge lags behind the high temperature slow shock. The kinematical estimations of the 
surge (Fig. 3) also reveals the acceleration of the surge material and its subsequent free fall 
due to gravity thereafter, which also matches well the simulation results. However, some 
parts of the comparison between the numerical and observational data are approximate 



- 15 - 



and only qualitative. This mismatch may result from simplified profiles of the parameters 
in our numerical model. Moreover, the real surge was excited in more complex plasma and 
magnetic field conditions at the boundary of AR11271. 

It is reported in literature that if the height of the energy release site is less (greater) 
than a critical height then the jet will b e driven by crest sh ocks (directly by the pressure 



gradient force generated by explosion) ( jShibata et al. 



19821 ). However, this minimum 



critical height in the solar atmosphere depends upon the strength of pressure enhancement 
in time at the energy release site. For the larger pressure enhancement, this height will 
be lower and vice versa. Under the conditions when the energy release site is higher in 
the solar atmosp here, the sur^e material may be accelerated by t he J xB force as a resul t 



of reconnection (lYokoyama fc Shibata 



1995, 



1990) as observed by 



Nis hizuka et al. 



torn . 



Therefore, the surge dynamics may not be so simple like a thermal pulse driven jet, and 
magnetic forces may also play some role in its formation. If the reconnection and thus 
the energy release site occurs in the lower chromosphere, any jxB force driven fiow may 
generate slow mode MHD wave /sho ck which eventually accelerates the jet in the upper 



chromosphere and transition region ( IShibata et al. 



19821 ). We consider the energy release 



site and the triggering by the thermal pulse at ^=1.75 Mm above the photosphere, and its 
amplitude Ap=80 times large r than the ambient R as pressure. In our case, is beyond 

( 119821 ). We aim to pe rform some pa r ametr ic 



the parameters considered by 



Shibata et al. 



studies in order to compare numerical results with the findings of 



Shibata et al. 



(Il982b, 



which will be devoted to our potential future studies. However, it is noteworthy tha t the 



numerical model we devised, is an extension of the ID model of Shibata et al. 



(119821). We 



adopted the temperature profile and a curved magnetic field in the MHD re gime, therefore 



our 2D model represents a natural extension of the hydrodynamic model of 



Shibata et al. 



(119821 ) in simulating the cool jet. 
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Table 1. The Summary of Input Free Parameters 



Half Threshold level Min filling Output Minimum Max traced Min curvature 
Width in standard factor resolution length structures radius 

flux deviation 



To trace the surge 2 pixels 1.0 0.35 5.0 0.01 10,000 25 

To estimate the plasma properties 10 pixels 1.0 0.15 5.0 0.04 10,000 60 
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Fig. 1. — SDO/AIA 193 A (right, yellow-brown map for 1.0 MK plasma) and 304 A (left, red- 
temperature map for 0.1 MK plasma) EUV images overlaid by HMI magnetic field contours of 
maximum level ±800 G. Yellow (blue) contours show positive (negative) polarity of the magnetic 
field. The surge origin is accompanied by heating above its foot-point. The observed surge is 
clearly evident in He II 304 A filter that is sensitive to the low temperature plasma of 0.1 MK. 
Some heated part of the surge especially near its base up-to the height of 20 Mm, is also evident in 
the coronal filters (e.g., 193 A ) that are sensitive to the plasma temperature around 1.0 MK. The 
surge is originated from the western boundary of the positive polarity situated in the North- West 
direction of the big sunspot with negative polarity as represented by the blue circular contour at 
X = 700", Y = 200". 
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Fig. 2. — SDO/AIA 304 A time sequence of a surge propagation higher into the corona from 
the boundary of an active region NOAA AR11271. The surge fades graduahy with time, and the 
material tries to settle down due to gravity. Horizontal and vertical dotted lines on the snapshot 
of 09:21 UT show the two paths along which the Distance-Time measurements have been made. 
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Fig. 3. — The height-time plots : (i) Along the solar surge using He II 304 A data (top) where the 
distance is measured along the slit chosen along the surge; (ii) Across the surge at a height of f»15 
Mm from its base (bottom) where the distance is measured along the slit chosen across the surge. 
The positions of both the slits are shown over the 09:21 304 A snapshot in Fig. 2. 



Fig. 4. — The Emission Measure (EM) map (top-left) and temperature map (top-right panel) : 
Evidence of the multi-temperature plasma forming the surge. The lower part of the surge consists 
of cool plasma, while the leading edge exhibits coronal temperatures. Bottom panel shows the 
automated tracing of the plasma threads as visible at various temperatures along the o pen field 



lines of the surge as shown by arrow using the methodology of 



Aschwandenetal 



mm . Colour 



codes for the emitting plasma in different AIA channels are shown at the top of this snapshot. 
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Fig. 5. — The AIA 171 A image (top-panel) displays the sub-field containing the surge and its 
overlying diffuse corona. The image axes are scaled as a fraction of the solar radius. This panel 
shows the 171 A image in which the box shows the area along the surge and its overlying diffuse 
corona along which the physical parameters are evaluated. The estimated plasma properties, e.g., 
the integrated density, temperature along the surge (dashed-line) during its maximum rise are 
displayed in two consecutives middle-panels. 
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Fig. 6. — Co-aligned SDO/HMI time sequence of the photospheric magnetic field near the origin 
of the surge. The blue surge is shown as intensity contour on HMI snapshot of 09:21 UT. The blue 
loop arches on both the sides of the surge may cause possible reconnection at its base to heat and 
subsequently launch a thermal pulse for its acceleration. 
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Fig. 8. — The grid system with block boundaries represented by sohd hnes as adopted in the sim- 
ulation. The initial pressure perturbation of the system is displayed by the horizontally elongated 
patch at a location (OMm, 1.75Mm). 
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Fig. 9. — Temporal snapshots of a simulated surge. The temperature (colour maps) profiles at 
t = 200 s, t = 400 s, t = 600 s, t = 800 s, t 1000 s, and t 1200 s (from top to bottom). 
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Fig. 10. — Temporal snapshots of a simulated surge. The density (colour maps) profiles at t = 200 
s, t = 400 s, t = 600 s, t = 800 s, t = 1000 s, and t = 1200 s (from top to bottom). Density is 
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